Understanding critical behavior in the framework of the 
extended equilibrium fluctuation theorem 
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Abstract. Recently, we have derived a fluctuation theorem for systems in thermodynamic equi- 
librium compatible with anomalous response functions, e.g. the existence of states with negative 
heat capacities C < 0. In this work, we show that the present approach of the fluctuation theory 
introduces new insights in the understanding of critical phenomena. Specifically, the new theorem 
predicts that the environmental influence can radically affect critical behavior of systems, e.g. to 
provoke a suppression of the divergence of correlation length £ and some of its associated phenom- 
ena as spontaneous symmetry breaking. Our analysis reveals that while response functions and 
state equations are intrinsic properties for a given system, critical behaviors are always relative 
phenomena, that is, their existence crucially depend on the underlying environmental influence. 
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1. Introduction 



Critical phenomena is the collective name associated with the physics of critical points pQ, which 
is originated as a consequence of the divergence of the correlation length £. The understanding of 
critical phenomena has motivated the development of new mathematical tools for the study of phase 
transitions; such as the renormalization group theory, which has a significant impact in condensed 
matter and high energy physics [2 . In this work, we present simple theoretical arguments and 
simulations that reveal new aspects in the understanding of critical phenomena. We show that the 
environmental influence can radically affect the critical behavior of systems, as example, to provoke 
a suppression of the divergence of correlation length £ and some of its associated phenomena as 
spontaneous symmetry breaking. 

Our arguments are based on a framework of equilibrium fluctuation theory recently proposed 
[21 IH 111 [7] . The main contributions of this approach are the derivation of a set of fluctuation 
theorems compatible with anomalous response functions [3j [4j 13 [6] as well as uncertainty relations 
involving thermodynamic quantities [7]. For the sake of self-consistence of this paper, we firstly 
present a brief review of the so-called fundamental equilibrium fluctuation theorem and some of 
its most immediate consequences [6]. After, this theorem will be considered to analyze the role of 
the environmental influence on critical phenomena. The results of this analysis will be illustrated 
studying the critical behavior 2D Ising model. Finally, we present some concluding remarks. 



2. Extended equilibrium fluctuation theorem 

The fundamental equilibrium fluctuation theorem (EFT): 

K = C + 1ZD (1) 

describes the relation between the system response functions and its fluctuating behavior [6]. The 
quantity TZ represents the response matrix: 

_ ( d p H dp (fiX) \ , , 

K ~ \ d Y U pd Y X ) ' [Z> 

which characterizes the system response to the environmental influence; whereas Y = 
(p, —E, —H, . . .J denotes a generalized force (pressure, electric field, magnetic field, etc.) and X = 

(v, P, M, . . .J the corresponding generalized displacement (volume, polarization, magnetization, 

etc.). Besides, H. = U + YX and U are the system Enthalpy and the internal energy, respectively; 
while f3 = 1/T0 denotes the inverse temperature and d x A, the partial derivative dA/dx. The 
quantity C is the self-correlation matrix that characterizes the system fluctuating behavior: 

/ (SQ 2 ) /3 (SQSX) \ 

\P (SXSQ) (3 2 (6X 2 ) J ' [6 > 
while the correlation matrix D: 

D= f (5P6Q) P(6P5X) \ 



(SYSQ) 13 (SYSX) 

describes the existence of environmental feedback effects among the system macroscopic observables 
(U, X) and the environmental control variables (/3, Y) due to the underlying thermodynamic 



| Boltzmann constant fcg has been set as the unity. 
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interaction. As expected, the amount of heat exchanged between the system and the environment 
at the equilibrium SQ = SU + YSX = TSS obeys the condition (SQ) = 0. 

EFT (JlJ represents a suitable extension of the usual fluctuation theorem [8]: 

TZ = C (5) 

derived from the Boltzmann-Gibbs (BG) distributions: 

d PBG (U,X\f3,Y) = eXP[ ~z%Y) YX)] m X ^ dUdX - ( 6 ) 

Here, ft(U,X) is the system density of states, while Z(f3,Y) is the partition function. Eq.([5]), as 
example, contains the fluctuation relations: 

C p = (3 2 (6Q 2 ) ,VK T = p (SV 2 ) , xt = (3 (SM 2 ) , (7) 

where C p = T(dS/dT) p is the isobaric heat capacity, Kt — —V^ 1 (dV/dp) T is the isothermal 
compressibility and xt = [pM /dH) T is the isothermal magnetic susceptibility. Despite the wide 
application in statistical mechanics, BG distributions © have a restricted applicability due to 
the environmental inverse temperature /? and its generalized force Y are regarded as constant 
parameters. Such that restriction demands the consideration of a thermal contact with a bath 
with infinite heat capacity or a reservoir with an infinite number of particles, as other idealizations 
assumed in conventional applications (usually attributed to the natural environment). In a general 
equilibrium situation, the internal state of the system acting as "environment" is also perturbed 
by the underlying thermodynamic interaction. This is the origin of environmental feedback effects 
described by the correlation matrix D in the extended EFT ([1]). 

Conventional EFT ([5]) is only compatible with states with a positive definite response matrix 
1Z. Considering its particular fluctuation relations ([7|. this last requirement implies the positive 
character of the heat capacity C p , the isothermal compressibility Kt and the magnetic susceptibility 
Xt- Remarkably, thermodynamics also supports the existence of states with negative heat capacities 
and other anomalous response functions, that is, a non-positive definite response matrix 1Z. 
Mathematically, the existence of anomalous response is associated with the presence of states where 
the entropy S = log W is a non- concave function. Such a relationship comes from the possibility of 
expressing response matrix 1Z as follows: 

■R = -TBT l T T , (8) 

where H is the entropy Hessian: 

/ d 2 S dudxS \ 

dxduS d 2 x s J ' w 

and the transformation matrix T: 

Moreover, the thermodynamic quantities /3 and Y are also obtained from the entropy S via the 
expressions: 

13 = (dS/0U) x , (3Y = (dS/dX) v . (11) 

The states with anomalous response conform the so-called regions of ensemble inequivalence 
[21 31 O H] . Physically, the anomalous response is a consequence of the incidence of non- extensive 
effects [31 HI [51 [3] , as example, the presence of long-range interactions or the development of spatial 
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non-homogeneities during the occurrence of discontinuous phase transitions [91I10|[TT|. The presence 
of a long-range force as gravitation explains the existence of negative heat capacities in astrophysical 
systems [9| I10| fTTj. The development of interphases induces the existence of surface correlations, 
which are the origin of the negative heat capacities C < observed during the phase coexistence 
phenomenon [9] . Analogously, such surface correlations explain the anomalous states with negative 
susceptibility \t < observed in a magnetic system below the critical temperature, which owe their 
origin to the formation of magnetic domains [6]. 

Direct consequences of the conventional EFT ([5]) are known [8] : 

• The system fluctuating behavior is determined by its own response to the external conditions. 

• The stable states are those ones where the response matrix 1Z is also positive definite because 
of the self-correlation matrix C is always positive definite. 

• The simultaneous divergence of some components of the response matrix and their 
corresponding correlation functions at the critical points, TV J — > oo C lJ — > oo, which is 
related to the divergence of correlation length, £ — >• oo. 

The presence of the correlation matrix D u in the extended EFT (Q} introduces the following radical 
changes in the previous implications: 

• The system fluctuating behavior depends both on its response functions and the nature of the 
environmental influence. 

• The stable states are those ones where 1Z(I — D) is a positive definite matrix, with / being the 
identity matrix. States with anomalous response functions can be thermodynamically stable 
under certain environmental conditions. 

• The divergence of the correlations functions C tJ and the correlation length £ can be suppressed 
by the environmental influence despite the divergence of some components of the response 
matrix TV^ at the critical point. 

The first and the second implications, together, claim that both the system fluctuating behavior 
and its thermodynamic stability depend on the nature of the environmental influence. A trivial 
example is found in the case of the isolated system, where SQ = 0. Clearly, there is no direct 
relation between the heat interchange and the heat capacity C p for this particular situation, which 
gives evidences about the restricted applicability of the fluctuation relation C p — 1 (<5Q 2 ). For a 
better understanding, let us consider a particular system where the unique relevant observable is the 
internal energy U. In this case, the EFT (JTJ drops to the following energy-temperature fluctuation 
relation [3]: 



Eq.(|12p is compatible with the existence of states with negative heat capacities C < as long as 
the correlation function (Sf36U) obeyed the stability condition: 



The simplest scenario with non- vanishing correlated fluctuations, (6/35U) ^ 0, is where the system 
is put in thermal contact with other short-range interacting system with a finite heat capacity C e . 



C = f3 2 (SU 2 ) + C (SPSU) , (12) 

which constitutes an extension of the fluctuation relation C = (3 2 (SU 2 ) obtained from the Gibbs' 
canonical ensemble (a special case of BG distributions): 




(S(35U) > 1. 



(14) 
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In this case, the second system (environment) experiences a temperature fluctuation 5T = —5U/C e 
whenever the system absorbs or releases the amount of energy SU. Substituting this last expression 
into Eq. (fT2|) . one obtains the self-correlation function of the internal energy U: 
i " / ' 

Stc = ? ^ ' (15) 

Accordingly, the environmental heat capacity C e drives the system fluctuating behavior and imposes 
the standard microcanonical and canonical conditions in the asymptotic limits C e — > and C e — > oo, 
respectively. The stability condition (fT4|) leads to the so-called Thirring's constraint C e < C 
necessary for the stability of negative heat capacities in the framework of short-range interacting 
systems [ID] . 

3. Implications on critical phenomena 

Let us enter into the central discussion of this work: the implication of the extended EFT (p} on 
critical phenomena. Our interest is to analyze the thermodynamic behavior of a magnetic system 
when the same one is immersed into an environment with constant inverse temperature /?; at the 
same time, it is put under the influence of an external magnetic field, H, that experiences a non- 
vanishing magnetic feedback effect (SHSM). This kind of environmental influence naturally appears 
when the magnetic field H s associated with the system magnetization M affects the source of the 
external magnetic field H and induces the existence of non- vanishing correlations (SHSM). For the 
sake of simplicity, it has been assumed that the external magnetic field H is weakly perturbed by the 
magnetization fluctuations, so that, this effect can be considered within the linear approximation 
H = Hq — XSM/N; where N is the system size, Hq = (H), the expectation value of the external 
magnetic field H, and A, an effective coupling constant that characterizes the system-environment 
magnetic interaction. In the case of A = 0, it considers the equilibrium situation described by BG 
distributions (J6J. Conversely, a value A ^ only ensures the constancy of the external magnetic 
field H in average sense. The application of the extended EFT ([JJ to this type equilibrium situation 
leads to the following relation for the susceptibility xt- 

f3 XT = (3 2 (SM 2 ) - (3 2 X t (SHSM) , (16) 

and therefore: 

Admitting an extensive growth of the susceptibility \t with N as xt oc N, it is possible to 
verify that the self-correlation of the external magnetic field behaves as (5H 2 ) oc 1/N. From 
the thermodynamic viewpoint, A ^ does not differ from A = when N — > oo. However, the 
positive character of the self-correlation function (<5M 2 ) now leads to the following condition of 
thermodynamic stability: 

A + N/xt > 0. (18) 

According to Eq. (ll7[) . the self-correlation function of magnetization (SM 2 ) remains finite when 
Xt — ► oo as long as the coupling constant A takes a positive value. Consequently, there is 
no divergence of the correlation length £ for this special equilibrium situation. Moreover, the 
consideration of a coupling constant with A > enables the access to diamagnetic states xt < 
that exist below the critical temperature, which are schematically represented in Fig. [TJ 
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Figure 1. Panel (a) Schematic representation of the regions of thermodynamic stability of a 
magnetic system within the BG distributions l(6]) in the M — U diagram (magnetization-internal 
energy). White-region: stable states with \T > an d M ■ H > 0; Gray region: metastable states 
with XT > and M ■ H < Oj Light-gray region: unstable states with \T < and M ■ H < 0. 
The critical point C with = /3 C and H = is a state of marginal stability for BG distributions. 
Panel (b) Behavior of some typical isotherms in the H — M diagram. 



Let us suppose that the magnetic system is initially set at a state O in the line with zero 
magnetization at Fig. [TJa, where j3 < f3 c and Ho = 0. Considering that the environmental inverse 
temperature (3 is increased with and without the magnetic feedback effect (SHSM), a relevant 
question emerges: What is the system behavior under these two external influences? The states with 
zero magnetization, M = 0, inside the region with xt < are unstable within the BG distribution 

In this case, the self-correlation function (6M 2 ) diverges for A = when the system approaches 
the critical point C (where xt — > oo). For j3 > f3 c , the system is forced to move along the stable 
symmetric curves a or a' with non-vanishing magnetization, where it undergoes a spontaneous 
symmetry breaking. Conversely, diamagnetic states \t < turn stable choosing a positive value 
of the coupling constant A. In this latter situation, the self-correlation function does not 

diverge when the system approaches the critical point C. For (3 > f3 c , the system remains on the 
line with zero magnetization M = inside the region with \t < 0, which implies a suppression 
of the spontaneous symmetry breaking. Thus, the system critical behavior can be modified by the 
environmental influence. Noteworthy that the incidence of non-vanishing feedback effect (SHSM) 
can also force the divergence of the self-correlation function ((5M 2 ) at an arbitrary state P inside 
the region with negative magnetic susceptibility. The occurrence of this phenomenon demands 




Figure 2. Inverse temperature dependence of the thermodynamic behavior of the 2D Ising model 
with L = 25 and Hq = with (open circles) and without (black circles) a magnetic feedback effect 
(5H5M) obtained from the Monte Carlo simulations with n = 10 5 iterations per point. 



a coupling constant A = —N/xt (P) > 0. Thus, the simultaneous divergence of the correlation 
length £ and the response function \t is only admissible if A = 0, that is, in the framework of BG 
distributions ©. 

The previous analysis can be verified through appropriate Monte Carlo (MC) simulations. 
For this purpose, we consider the Ising model on a square lattice L x L with periodic boundary 
conditions, whose internal energy U = — Y)/jj\ &iO~j and magnetization M = <7j, where cr, = ±1 
and (ij) denotes the nearest neighbors sites. As the case of conventional Metropolis importance 
sample [13] based on the BG distributions ©, the present situation with (5H5M) ^ can be 
implemented through the following acceptance probability p for a MC step from the state (U, M) 
to (77 + AU,M + AM): 

p = min {exp [-p (At/ + HAM)] , 1} , (19) 

where the magnetic field H has an implicit dependence on the magnetization fluctuations. The 
magnetization dispersion = (<5A/ 2 ) /N can be employed to obtain the inverse susceptibility per 
particle as x^ 1 = N\t 1 — (l — ^P^m) I ' P a m- Along the Monte Carlo simulations with {5H5M) ^ 
0, the value of the coupling constant A is redefined for each calculated point using the optimal value 

X op t = \J 1 + (x^ 1 ) 2 — Xt 1 ' wmcn minimizes the total dispersion a 2 = (<5M 2 ) /N + N (SH 2 ) and 
satisfies the stability condition (|18l) . 

A comparative study of Monte Carlo simulations with and without a magnetic feedback effect 
(SHSM) is shown in FIG12] According to results that we show in FIGl2]c, the magnetization 
dispersion diverges at the critical point (3 C {L = 25) ~ 0.41 in Monte Carlo simulations with 
coupling constant A = 0. Conversely, the same quantity remains finite when A > 0, which 
indicates a suppression of the long-range order associated with the divergence of correlation length 
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£ — » oo. Simultaneously, there exist a suppression of the spontaneous symmetric breaking, which is 
manifested in FIG|2]a and FIGjSJb as a persistence of diamagnetic states, \t < 0, with (M) = for 
(3 > f3 c . Remarkably, the absolute values \xt\ of these negative susceptibilities are very large. The 
suppression of the spontaneous symmetric breaking is also manifested as a bifurcation of caloric 
curves and the inverse isothermal susceptibility x^ 1 at the critical point C shown in FIG Ob and 
FIG(2]d, a fact that gives evidences to the different thermodynamic behavior of this system inside 
the region with negative susceptibilities. 

4. Concluding remarks 

As already evidenced, the extended EFT ([T]) reveals some new insights in the understanding 
on statistical physics. This theorem clearly distinguishes two different type of thermodynamic 
properties for a real system: the intrinsic properties and the relative ones. Examples of intrinsic 
properties are the response junctions ([2]), the state equations as well as any other information 
derived from the thermodynamics of the isolated system; in particular, from the knowledge of 
the entropy S. Conversely, the relative thermodynamic properties depend on the environmental 
influence acting in a concrete equilibrium situation, as example, any information concerning to 
the system stability and its fluctuating behavior. The divergence of correlation length £ and its 
associated critical behaviors (spontaneous symmetry breaking, critical opalescence, etc.) take place 
at states with a marginal stability (see in Fig[TJa). Since stability conditions, e.g.: Eq. (|18[) . crucially 
depend on the environmental influence, any critical behavior is always a relative phenomenon. 
Hence, the divergence of any response function component 7U J — > oo is not an authentic critical 
behavior, but an intrinsic phenomenon that has nothing to do with the divergence of correlation 
length £. The above remarks evidence that the consideration of BG distributions ([6]) can potentially 
lead to some incorrect predictions about the system thermodynamical behavior; overall, when one 
assumes these equilibrium distribution functions outside the context of its traditional applications. 
This could be the case of collective phenomena involving systems with long-range interactions as 
non-screened plasmas and astrophysical systems [TT], or, the small and mesoscopic systems as the 
case of nuclear multifragmentation 12] , where the existence of states with anomalous response is 
almost a rule rather than an exception. 

Before end this section, we would like comment three possible application frameworks of our 
approach, which could deserve a special attention in future works. Previously, this framework 
of fluctuation theory have been successfully applied to overcome slow sampling problems in 
Monte Carlo simulations associated with the incidence of discontinuous phase transitions O [15] . 
Analogously, the suppression of divergence of correlation length £ at the critical point could be 
also employed to design new Monte Carlo algorithms to deal with difficulties associated with 
continuous phase transitions. The second application framework are the quantum theories at finite 
temperature [16] . which exploit the analogy between the evolution operator T = e ~ ltH / h and the 
statistical operator ui — e~^ H ^ kB . In high energy physics, these formulations have been applied 
to study processes that took place at the early universe as a whole, as example, spontaneous 
symmetry breaking of electro-weak interactions. Since these processes are special cases of critical 
phenomena, the application of present framework of fluctuation theory could reveal new insights in 
their understanding. Finally, the analogy between electric and magnetic systems strongly suggests 
the existence of anomalous dia- electric states (electric counterpart of diamagnetic states) below 
critical temperature of a ferro-electric system, which should not be detected using experimental 
techniques based on Boltzmann-Gibbs distributions ([6]). We think that the ideas discussed in this 
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work could inspire new experimental techniques to clarify the existence (or nonexistence) of these 
anomalous states. As discussed elsewhere [17], the presence of a hypothetic dia-electric medium 
favors the formation of Cooper-pairs, a mechanism that triggers the development of high-temperature 
superconductivity. 
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